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ABSTRACT 


The equation defining the energy radiated per unit solid angle due to Cerenkov 
radiation is analyzed in detail, including the effects of all equation variables for a 
hypothetical electron accelerator experiment. Specifically, various finite frequency 
intervals are used in an effort to determine the optimum means of determining the 
details of a charge bunch in a high energy electron accelerator. In particular, it is 
shown how narrowband measurements as a function of angle may vield information on 
both the beam path length and the bunch charge parameters. As an aid to the 
analysis, an interactive Fortran program 1s presented that allows for any specific 


experimental parameters, with options for various output types as desired. 


THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may 
not have been exercised for all cases of interest. While every effort has been made, 
within the time available, to ensure that the programs are free of computational and 
logic errors, they cannot be considered validated. Any application of these programs 


without additional verification is at the risk of the user. 
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I. INTRODUCTION 


A. HISTORY 

The existence of the phenomenon now known as Cerenkov radiation was 
observed as early as 1910, most notably by Madam Curie. However, during the 
ensuing vears, there was other work being completed by those familiar with Madam 
Curie’s observation that masked a detailed study of the phenomenon. Jelley [Ref. 1] 
explains their reasons in greater detail. In 1926, Mallet made the first deliberate 
attempt to study and explain “the very faint emission of a bluish-white light from 
transparent substances”. Experimentation concerning this phenomenon continued 
through the 1920’s and 1930’s by Cerenkov and Mallet, but a viable theory explaining 
the process would not be proposed until 1937 by Frank and Tamm. Their theory, for 
which they were consequently awarded the 1958 Nobel Prize in Physics, was found to 
be in excellent agreement with the experimental results obtained by Cerenkov. The 
advent of more sensitive light detectors and other experimental equipment accounted 
for more in depth studies of Cerenkov radiation in the 1940’s and 1950’s. [Ref. 1: pp. 
1-8] 

Since then, numerous experiments have been conducted and many papers and 
academic theses published that have helped to more clearly understand and explain the 
Cerenkov radiation phenomenon. Some of the more recent investigations include; 
periodic electron bunches of finite emission lengths, time developement of bunch 


charges and Cerenkov radiation in the x-ray region, (see Bibliography). 


B. BACKGROUND 

Work on this thesis was motivated by two factors. First, using previously derived 
equations and expressions describing the energy radiated per unit solid angle due to 
Cerenkov radiation, [Ref. 2] a fortran program was written to run on the IBM 3033 
mainframe. The program, (see Appendix), was written interactively to allow all 
parameters to be changed in order to tailor the output to fit any specific experimental 
setup. 

Second, using the aforementioned Fortran program, a series of outputs were 
generated to be used in comparing the theoretical results of previously derived 


expressions with those of actual experimentation. Since there are literally an infinite 


number of combinations of various parameters, only a few were chosen to be included 
in this work. The parameters chosen are consistent with those expected while 


conducting experiments on any high energy electron accelerator. 
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It. THEORY 


Aww) CBRENKOV EFFECT 

Cerenkov radiation can most easily be explained by describing the interaction of 
a single electron within a transparent medium. Consider a non-relativistic electron. 
While traversing the meditim, the electron will tend to locally polarize the adjacent 
atoms in the medium, instantaneously creating a temporary dipole within the material. 
Therefore, as the electron moves through the medium, an electromagnetic pulse is 
generated. Ilowever, since there is complete symmetry, there will be no net field 


Beneratcd, Ste figure 2.1. 
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Figure 2. Polarized Atoms in a Dielectric. 


Now suppose an electron ts travelling at a speed comparable to or exceeding the 
speed of light for the medium. Although symmetry 1s preserved in the azimuthal plane, 
along the axis of motion a resultant dipole field 1s generated. These fields are set up at 
each element along the electron’s track, radiating a brief electromagnetic pulse. 
Provided the electron is moving at a speed greater than the speed of light within the 
medium, the wavelets formed from all elemental positions on the track can be in phase 
and thus produce a resultant electromagnetic field. [Ref. I: pp 3-6] 

igure 2.2 depicts the relationship of the distance (BZ) travelled by an electron in 
time At versus the distance (B-1) covered by an emitted electromagnetic pulse during 


the same time. The distance (BZ) travelled by the electron is given by: 


i 


BZ Sacre (eqne2. |) 


Where f} is the ratio of the electron speed to the speed of light in a vacuum Cy The 


distance (2.1) travelled by the electromagnetic wave is given by: 

BA = cAt, (eqn 2.2) 
where c is the speed of light for the medium. For a medium with an index of refraction 
n, c=cC)/n. From equations 2.f and 2.2, the “cerenkov relation” is obtained. 


[Ref. I: pp 6-8] 


cosO. = 1/Bn. (eqn 2.3) 


VAt 





Pigyre 2.2 Illustration of Cerenkov Radiation. 


Due to the finite path length at which the radiation is detected, the Cerenkov 
cone angle is normally shifted away from 9 ~ As the frequency at which the emitted 
radiation 1s measured or the distance from the source is increased, the radiation cone 
angle approaches OF As this occurs, an increasing fraction of the total Cerenkov 


radiation is found at @.. (Ref. 2] 


B. MATHEMATICAL INTERPRETATION 
1. Energy Equation 
As previously discussed, for a finite path length, diffracted Cerenkov radiation 
effects will be seen at angles other than the Cerenkov angle. Previous work on this 
phenomenon [Ref. 2] has resulted in an expression for the emitted Cerenkov radiation. 
The energy radiated per unit solid angle within the frequency range dv by a single 


bunch charge q travelling a finite distance L 1s: 
E(v,k)dv = QR?dv, (eqn 2.4) 


where Q is a constant whose magnitude is dependent on the index of refraction and the 


total charge in Coulombs. 
Q = pWcq*/8r". (eqn 2.5) 
The radiation-function R is given by: 
R = 2nysinOI(u)F(4), (eqn 2.6) 


where @ is the angle measured from the line of travel of the bunch charge to the 
direction of propagation of the emitted radiation, I(u) is the diffraction function, and 
F(k) is a form factor. For convenience, the length of travel of the bunch charge is 
measured in units of the wavelength of the emitted radiation within the medium. The 
dimensionless beam length parameter nN serves this purpose and 1s defined as the ratio 
of the length L from the source to the wavelength A of the emitted radiation. 


n= L/A. (eqn 2.7) 
The diffraction function is defined as: 
I(u) = sin(u)/u, (eqn 2.8) 
where u is dependent on both the angle @ and the beam length parameter Nn. 


u = mn{(1/nB) - cos8}. (eqn 2.9) 
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Equations 2.5 through 2.9 can then be used to redefine equation 2.4 in terms of 
experimental constants, (ie. n, L, f etc.) and the single remaining variable of 


frequency. [Ref. 2] 
E(v,k)dv = QG*F?(k)sin?(Av}dv. (eqn 2.10) 


For a given experimental setup, A and G are functions of the angle only and are 


defined as follows: 


A = {Lm/c}e{(1/nB)-cos@}, and (eqn 2.11) 


G = {2sin9}/{(1/nB)-cos®}. (eqn, 2aien 


The radiation patterns arising from equation 2.4 can be thought of as an 
oscillating sin°x function modulated by an envelope. If F(X) is neglected, then G* acts 
as the modulating envelope. The actual form of the envelope will vary with nB, but 
will be constant for a given experimental setup. Figure 2.3 shows a plot of G? as an 
illustrative example for values of nB two percent above and below threshold, 
(nB = 1.0+0.02). 

Finally, the dimensionless form factor F(x) is related to the Fourier transform 
of the bunch charge distribution. The Fourer components of any bunch charge are 


defined as: 

p(k) = [p(r)ea?r, (eqn 2.13) 
where p(r) defines the charge distribution of the bunch charge. Once the Fourier 
components are determined from equation 2.13, the form factor is found by the 
relation: 

p(k) = qF(&). (eqn 2.14) 
Here we consider only a line charge with distribution p(z) which is travelling in the z 


direction. Thus, r-+z in equation 2.13 and k-+k, in equation 2.14. In this work, only 
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Figure 2.3 Illustrative Example of the Cerenkov Radiation Envelope. 


the trapezoidal bunch charge distribution is used although many others are available 
and easily used within the Fortran program. [Ref. 3] 
2. Trapezoidal Form Factor 
The first step in obtaining any type of form factor is to determine the assumed 
geometry of the emitted bunch charge. Figure 2.4 depicts the geometry for the 
trapezoidal bunch charge. The functions defining the positive half of the trapezoid are 


given by: 


f(x) = A, and (eqn 2.15) 


f(x) = {A/(d-b)} e{d-x}, (eqn 2.16) 


Where f, defines the trapezoid from 0 to b and f, defines the trapezoid from b to d. 
The amplitude of the bunch charge 1s denoted as A with units of coulombs per meter. 


The values b and d have units of meters and are easily obtained from the pulse length. 





figure 2.4 Trapezoidal Bunch Charge. 


The fourier components of this trapezoid are defined by equation 2.13 where 
p(k.) 1s the sum of the two functions f(x) andf,(x), (ie. the charge distribution of a 
single bunch charge), and k, is the component of the wave vector &. The details of 


solving equation 2.13, though tedious, are rather straightforward. (Ref. 3] 


Equation 2.17 is the solution to equation 2.13 and provides the Fourier components 


in the case of the trapezoid. 
p(k,) = { (2A)(d-b)/k*, }#{ cos(k,b) - cos(k,d) }. (eqn 2.17) 


The value of A is found by equating q, the total bunch charge, to the integral 


of the charge distribution. 
a OO 
ee pizjdz = A(b+d). (eqn 2.18) 


By substituting the results of equations 2.17 and 2.18 into equation 2.14. the form 


factor for the trapezoidal charge distribution is: 
F(k,) = { 2q/(d*-b*) }@{1/k*} #{ cos(k,b)-cos(k,d) }. (eqn 2.19) 


If the medium is assumed to be non-dispersive, then @=ck, and remembering 
that k= kcos®, the form factor can be written in terms of the angular frequency 
2nv=@. After substituting, the explicit relation for the trapezoidal form factor is 
given by equation 2.20 as follows: 


F(k,) = { c?/2n*(d?-b*) }@{ cos(Bv)-cos(Dv) }e{ I/v? }, (eqn 2.20) 
where B and D are given by equations 2.21 and 2.22 respectively. 


B = 2nbcos@/c, and (feqme2. 21) 


D = 2ndcos@/c. (eqn 2.22) 


It is clearly seen that once the parameters for the bunch charge are defined, 
the form factor is highly dependent on the inverse square of the frequency v. The 
expression for the energy radiated per unit solid angle within the frequency range dv 
Written in terms of frequency will then be equation 2.4 with the form factor as given in 
equation 2.20. [Ref. 4] 


i 


Edv = [ QG?sin?(Av) ] iv 
*[ (c7/277(d?-b7)} {cos(Bv)-cos(Dv)} {1/7} }2dv. 
The radiated energy is a function of frequency through both the sin*(Av) term and the 
form factor. 


(eqne 2-23) 


II. RESULTS AND ANALYSIS 


The Fortran program written for this work has various output capabilities. The 
first type of output contains data equating the explicit value of E (eqn.2.4), in joule- 
seconds, to a specific frequency at a particular angle 8. Another output type provides 
data for graphing the radiation envelope (eqn.2.3) or the form factor (eqn. 2.20) used in 
this thesis. These graphs are useful in understanding the complex nature of the 
variables involved in solving equation 2.4. 

The most useful output type is that which solves, (1.e. integrates), equation 2.4 
over a predetermined frequency range for all desired angles from the beam path. A 
careful study of these graphical outputs will provide insight into the particulars 
necessary for any experimental setup so that the most accurate and informative results 
can be obtained. 

As previously discussed, the output obtained during this work are based on a set 
of parameters that uniquely define a specific electron accelerator experimental setup. 
The following input parameters (with their Fortran variable names), are required: 

1. Accelerator Beam Energy (EBEAM). 
Total Periodic Electron Bunch Charge (CUE). 
Index of Refraction (IND). 
Distance from Source to Detector (LENGTH). 
Cerenkov Radiation Frequency Interval (MIDNU,ENU). 
Incremental Change in Frequency Interval (DNU). 
Angle of Interest, as Measured from the Charge Path (DITHETA). 
Bunch Charge Form Factor (FORFAC). 
Trapezoidal Bunch Parameters (BEE,DEE). 


The specific values for each variable used to obtain the outputs in this work are 


to 


. he ele he Bie 


as follows: 

1. Beam Energy - 50 MEV. 

2. Bunch Charge - .001 Coulombs. 
3. Index of Refraction - 1.000268 (air). 
| 


Source to Detector Distance - 100 Meters. 


Ca 


Frequency Interval - As noted. 


es 


Frequency Increment - 0.5 MHz. 


7. Solid Angle - As noted. 
8. Form Factor - Unity and Trapezoid. 
9. Trapezoidal Bunch Parameters. 
a. Top-50 nanoseconds. 
b. Base-60 nanoseconds. 
For the given index of refraction and beam energy, the “Cerenkov Relation” (eqn.2.3) 
gives the critical Cerenkov angle as 1.9°. 

Three basic assumptions were made in the course of this work. First, it is 
assumed that the shape of the trapezoidal bunch charge remains unchanged as it 
travels through the medium. This is reasonable since the distances involved are 
relatively short with respect to the speed at which the bunch charge is travelling. 
Second, the permeability of free space wt is used vice that of the actual medium of air. 

The final assumption is that the speed of the electron bunch remains unchanged 
throughout it’s travel. Although there is some Bremstrahlung radiation emitted, it has 
a negligeable effect on the speed of the bunch charge and therefore B is assumed to 
remain constant. 

Figure 3.1 depicts the energy radiated per unit frequency at angles of 30, 45 and 
60 degrees using a form factor of unity. As would be expected from equation 2.23, a 
sin*x function results since Q (eqn.2.5), A(eqn.2.11) and G(eqn.2.12) are all constant at 
the given angle 9, the only variable being the frequency v. The difference in 
amplitudes is due to the variation of G with @. The difference in the periods is due to 
the function A, which, as the argument of the sin term, effects the periodicity. 

Once a form factor other than unity is used, the output pattern is changed 
significantly. Equation 2.23 shows the frequency dependence of the trapezoidal form 
factor. Figure 3.2 is the form factor obtained at 0 degrees over the frequency range 
10-100 MHz. The damping of the form factor waveform is due to the inverse square 
relationship of F(X) with frequency. Figures 3.3 and 3.4 depict the square of the form 
factor as it is used in equation 2.4. Different frequency intervals are used but each 
waveform is computed at the same angle, 9=0°. As with Figure 3.2, Figures 3.3 and 
3.4 each dampen with increasing frequency. The increased periodicity of Figure 3.4 is 
due primarily to the wider frequency range. At angles other than 9@=0° the form 
factor will be similar in shape except scaled by changes in B and D due to the cos@ 
factor. [Ref. 4] 
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Figure 3.2 Trapezoidal Form Factor, (8 = 0°). 
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Figure 3.3 Square of the Trapezoidal Form Factor, 10-1l00MIIz., (8 = 0°). 
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Figure 3.4 Square of the Trapezoidal Form Factor, 200-800MIIz., (0 = 0°). 


Figure 3.5 is a plot similar to that of Figure 3.1 except that the unity form factor 
has been replaced with that of the trapezoid. As expected, the use of a specific form 
factor significantly effects the output. The output of interest will be the integration of 
all outputs similar to Figure 3.5 at each angle of interest. Before analyzing the 
integrated output with the trapezoidal form factor, it is useful to look at Figures 3.6 
through 3.8 which depict the integration of equation 2.23 while maintaining unity form 
factor. As noted earlier, an energy peak at the Cerenkov angle of 1.9° is expected. 
Although the peak appears to be centered at about 10° in Figure 3.6, as the frequency 
interval increases from 10-100 MHz. to 10-10,000 MHz. the same curve shifts closer to 
the Cerenkov angle. Concurrently, the graphs become smoother as more data points 
are used at each interval @. The aberration found in Figure 3.8 may be an artifact of 
the Fortran program. The extraordinarily large magnitude of Figure 3.8 is due to the 
existence of a pole in variable G (eqn.2.12) of equation 2.10 as @ approaches the 
critical Cerenkov angle. However, since the variable A is zero at the Cerenkov Angle, 
equation 2.10 remains finite. [Ref. 2] 

Figure 3.9 shows the output of equation 2.23 integrated over a relatively 
broadband frequency range of 100 to 1,000 MHz. Figure 3.10 is the identical output 
plotted over a narrower range of 8 to enhance the output detail. In each case, the 
detail of the output is fuzzy at best, making it difficult to conduct any meaningful 
analysis. In an attempt to force more detail in the output, the frequency range over 
which equation 2.23 1s integrated was narrowed. Figure 3.11 depicts the energy spectra 
integrated over 20 MHz. centered around 100 MHz. All other parameters are equal to 
those used in Figures 3.9 and 3.10. Although similar in general form , Figure 3.11 
shows more detail. Then, changing to a narrower frequency band of 2 MHz.(again 
centered around 100 MHz.), Figure 3.12 depicts a clearly defined output that is similar 
to a frequency modulated carrier wave. Figure 3.13 shows the same output for a 
narrower range of 8 to show more detail. 

The Cerenkov radiation envelope discussed earlier(Figure 2.3) has an impact on 
the outputs of Figures 3.12 and 3.13, but once nf is selected the envelope will effect 
only the amplitude of this and similar graphs. The null points of Figures 3.12 and 3.13 
can then be calculated by analyzing the components of the radiation function(eqn. 2.6). 
The expression for the trapezoidal bunch charge form factor(eqn.2.20) can be expressed 
as a product of a constant and two sin functions using the following trigonometric 


identity. 
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Figure 3.5 Energy Radiated per Unit Frequency, (8 = 30°, 45° and60°). 
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Figure 3.6 Energy per Unit Solid Angle, 10-100MEIz., (Unity Form Factor). 
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Figure 3.7 Energy per Unit Solid Angle, 10-l000MEIz., (Unity Form Factor). 
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Energy per Unit Solid Angle, 10-10,000MIHIz., (Unity Form Factor). 


cosa - cosy = -2sin{ (a/2)+(y/2) }esin{ (a/2)-(7/2) }. (eqn 3.1) 


After a great deal of algebraic manipulation, it can be shown that the form factor will 
be of the form: 


F(A) = const®.sin(€)®.sin(@), (eqn 3.2) 


where €= McosO(d+ b)v/c, and @ = ncosO(d-b)v/c. 

For the trapezoidal bunch charge parameters used in this work, c/(d-b) equates 
to a value of 200 MHz. and c/(d+b) to a value of 18.2 MHz. Then, using the mean 
frequency of 100 MHz. found in Figure 3.13, one need only to equate each of the sin 
function arguments to integer values of @ to mathematicallv determine the values of 0 
at which equation 3.2 will equal zero. Table 1 lists the values of 8 at which the form 


factor will force zeroes in the energy equation. 


TABLE 1 
VALUES OF 8 AT WHICH THE FORM FACTOR EQUALS ZERO. 


undefined. 
undefined. 
undefined. 
undefined. 


undefined. 





At the mean frequency of 100 MHz., sin @=0 cannot be determined since it 
requires solving for an angle @ whose cos is greater than unity. However, by 
comparing the values of 8 in column 1 of Table 1 to Figure 3.13, one sees a correlation 
of the null points to what can be considered the envelope of the output. It is clear 
then that the form factor acts as the modulator of the energy spectra of the Cerenkov 
radiation. 
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The “carrier wave” nulls can be determined in a similar manner by analyzing the 
diffraction function (eqn.2.8) found in equation 2.6.. Equating u to integer values of 
and solving for 8 will provide all angles 8 at which the energy output will be zero. 
Table 2 lists a sample of the results for which sinu equals zero. Again, comparing the 
values of 8 in Table 2 to Figure 3.13 , one sees that the calculated nulls equate to those 


of the “carrier wave’. 


Lee le, 
VALUES OF 8 AT WHICH SINU EQUALS ZERO. 





Figures 3.14 and 3.15 are graphs of the Cerenkov energy output over the same 
narrowband frequency interval of 2 MHz., but at the higher frequency range of 499 to 
501 MHz. Although some detail is lost compared to the outputs of Figures 3.12 and 
3.13, the same analysis can be conducted to determine the null points of both the 
“carrier wave” and the modulating form factor. Tables 3 and 4 list the values of @ for 
each null point similar to those found for Figures 3.12 and 3.13. The range of 8 was 
chosen so that a meaningful comparison could be made with the outputs of Figures 
3.14 and 3.15. 
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Figure 3.9 Energy per Unit Solid Angle, 100-1000MIIz., (Trapezoidal Form Factor). 
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Figure 3.10 Energy per Unit Solid Angle, 100-1000MHz., (Trapezoidal Form Factor). 
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Figure 3.11 Energy per Unit Solid Angle, 90-110MIIz., (Trapezoidal Form Factor). 
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Energy per Unit Solid Angle, 99-101 MH{z., (Trapezoidal Form Factor). 
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Figure 3.13 Energy per Unit Solid Angle, 99-101 MI1z., (Trapezoidal Form Factor). 


36 


TABLE 3 
VALUES OF @ AT WHICH THE FORM FACTOR EQUALS ZERO. 


undefined. 
undefined. 
undefined. 
undefined. 


undefined. 


TABLE 4 
VALUES OF 8 AT WHICH SINU EQUALS ZERO. 
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Figure 3.14 Energy per Unit Solid Angle, 499-501 MEfz., (Trapezoidal Form Factor). 
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Figure 3.15 Energy per Unit Solid Angle, 499-S501IMHz., (Trapezoidal Form Factor). 
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[V. DISCUSSION 


Figures 3.6 through 3.8 show that the energy radiated per unit solid angle from a 
charge distribution with unity form factor tends to be concentrated into an ever 
narrowing direction as the range of integration over frequency 1s increased. In the limit 
of all frequencies, the radiation is confined to the Cerenkov cone. The effect of 
including all frequencies is the same as allowing the beam path length to go to infinity, 
and the resulting radiation patterns are characteristic of that from a point charge. 

The Zeroes in the radiation patterns for trapezoidal charge distributions, Figures 
3.9 through 3.15, are tabulated in Tables | through 4. Since the tables and figures are 
calculated from the same theory, they are, of course, in agreement. Tables 2 and 4 
show the zeroes for the “carrier” oscillation which arises from the length of the electron 
beam path from source to detector. Tables 1 and 3 show the zeroes in the € and @ 
oscillations resulting from the length (d+b) and rise time (d-b) respectively of an 
individual charge bunch. Both Tables show zeroes dictated by the length of the bunch 
charge. Only two zeroes from the pulse rise time occur; both at the higher frequency 
of measurement. 

In the cases studied here, the bunch length (d+b) equals 16.5 meters and the 
pulse rise (d-b) equals 1.5 meters. At the measuring frequency of 100 MHz. the ratio 
of bunch length to wavelength is 5.5; large enough to resolve the details of the bunch 
size. Even at the higher frequency of 500 MHz. the ratio of pulse rise to wavelength is 
2.5 and consequently the rise and fall of the pulse have little effect on the radiated 
output. [Ref: 5] 
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V. CONCLUSIONS AND RECOMMENDATIONS 


After studying the figures depicting the radiated Cerenkov energy in this work, 
one could initially deduce that the use of narrowband detectors would be optimum for 
detecting radiation during any electron accelerator experiment. As shown and 
previously discussed, the information produced is more detailed when narrowband 
frequency ranges are used compared to wider band measurements. In particular, if one 
frequency band is considered as a function of angle, the calculated radiation energy has 
much detail as displayed in Figures 3.12 and 3.13 The structure depends on both the 
path length and bunch length, and, in principle, may yield information about these 
parameters. One must also consider the logistics of any experimental setup and any 
engineering problems involved to yield the maximum information. 

For the broadband output with a unity formfactor, the energy output has a peak 
at the critical Cerenkov angle as expected. This output is typical of a point charge. 
However, since the nominal electron accelerator bunch charge output has a finite 
distribution, there is little information to be extracted from this particular output. 

When the effect of the form factor is included, there exists a definite Cerenkov 
radiation peak perpendicular to the path the electron bunch 1s travelling. This 
phenomenon may be significant with respect to the optimum position at which the 
maximum energy produced by the radiation can be detected. 

If a narrowband detection system is realistically feasible, it can be used to 
determine all nulls that arise in the energy equation due to the diffraction function and 
the form factor. Although narrowband detection will produce the most detailed 
Output, its usefullness may be limited by the availability of equipment designed to 
discriminate the output within the desired frequency band limits. 

In order to obtain significant information about the details of the charge 
distribution within a bunch, it will be necessary to perform experiments at frequencies 
such that the ratio of mse time to frequency is at least 5. For real beams, where the 
rise times are of the order of a few nanoseconds, this means that narrowband 


measurements must be performed at microwave frequencies. 


4] 


To augment this research the following is recommended: 


|. Analyze outputs with form factors other than trapezoidal, (i.e. rounded, 
triangular etc.). 


2. Determine the optimum (higher) frequencv bands for various frequency 
intervals at which useful information can be extracted. 
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ArT ENDIX 
FORTRAN PROGRAMI 


PROGRAM TO COMPUTE THE ENERGY RADIATED PER UNIT SOLID ANGLE TRIOOOLO 
OVER A SPECIFIC FREQUENCY RANGE 


3633388 HE BEBE BEBE BBN BEEBE BEEBE BBEBEHEBE 
« ERG ARRAY SIZE MUST CORRESPOND TO THE NUMBER OF x 
x INCREMENTS OF THETA * 
JOB BBB EBB BBB BUBB RB ERB BBR BRB BBE BEEBE BHR SERB RB BHHHEE 


NgNANAANAANAAAN 


DIMENSION ERG( 1800 ) 
TRIOOO02O0 
TRIO0030 


ole) 


REAL ERG,COSANG ,»COSARG »COSGAR »>KKONST »DEE ,BEE 
10 REAL FORFAC ,EBEAM,LENGTH>KONST »>C ,»GEE, TOP ,BOT TRIO0040 
20 REAL ATHETA,RADFUN,PI,EREST »GAMMA,THETA,IND »MIDNU TRIOOO50 
30 REAL CONST ,MU,CUE,TERG,TTERG,DNU,BNU,DEGINC ,DATEND 

INTEGER IFORM,INCREM,DLOOP ,D2LOOP ,SELECT ,NUMBER 

CHARACTER*¥®1 RESPON,ANSH 


*%*x*XEBEAM=BEAM ENERGY (MEV) 5 FORFAC=FORM FACTOR (F(K)) 
*X**XETA=BEAM LENGTH/WAVELENGTH 3 BETA=VELOCITY/C( FOR MEDIUM) 
HHHGAMMA=RATIO OF BEAM ENERGY ELECTRON REST ENERGY 
x*x*DIFFUN=DIFRACTION FUNCTION (I(U)) § DTHETA=THETA( DEGREES ) 
RTHETA=THETAC RADIANS) 5 LENGTH=LENGTH FROM SOURCE( METERS) 
***XRADFUN=RADIATION FUNCTIONS MIDNU=BEGINNING FREQUENCY 
%**CUE=BUNCH CHARGE IN COULOMBSSEREST= ELECTRON REST ENERGY 
*#*KERGSARRAY FOR STORING VALUES OF RADIATED ENERGY 

REX FOR CUMULATIVE ANGLES THETA 
x**DNU=FREQUENCY INTERVAL $TERG=RADIATED ENERGY FOR A 

HEX SPECIFIC ANGLE THETA 


AAADANANAANAAAANAANAANAAAAADNA 


HHHASSIGNMENT OF CONSTANT S22 
CO=2.997925E8 
PI=3.14159 
EREST=.5117 
MU=1.2566E-6 


PRINT *,'THE INDEX OF REFRACTION (AIR) USED IS 1.000268' 
PRINT *,‘'DO YOU WANT A DIFFERENT VALUE?(Y/N)' 
READ (5,90 JANSH 
90 FORMAT(A1) 
IF (ANSW .EQ.‘'N' )THEN 
IND=1.000268 


GO TO 100 
ELSE 
PRINT *,'SELECT THE INDEX OF REFRACTION (REAL )' 
READ *,IND 
END IF 
100 CONTINUE 
C 
C DETERMINING THE SPEED OF LIGHT FOR THE MEDIUM 
C = CO/IND 
Cc 
PRINT *,'ENTER BUNCH CHARGE (REAL) IN COULOMBS' 
READ *,CUE 
260 CONTINUE 
c 
CCC COMPUTATION OF CONSTANT (CALLED @ IN THESIS) 
c 
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CONST =MUXC #( CUE #2. /(8. HC PI MHZ. )) 
CALL EXCMS ('CLRSCRN' ) 


PRINT *,°ENTER THE DISTANCE FROM THE SOURCE (IN METERS)’ 
PRINT *, ‘THE RECEIVER WILL BE PLACED’ 
READ *, LENGTH 
PRINT *, ‘SELECT THE FORM FACTOR DESIRED: ' 
PRINT *,°1=UNITY, 2=TRAPEZOID' 
READ *, IFORM 
IF (IFORM .EQ. 1LIJTHEN 
FORFAC=1.0 
GO TO 290 
ELSE 4 
PRINT *,‘ENTER THE TOP AND BASE VALUES (IN NANOSECONDS) FOR’ 
PRINT *,‘'FOR THE CURRENT TRAPEZOID FUNCTION' 
READ *,TOP,BOT 


Cc THE NEXT 2 LINES CONVERT PULSE TIME TO DISTANCE 
BEE C¥TOP/2.E9 
DEE CxBOT/2.E9 


END IF 
290 CONTINUE 
CALL EXCMS ( 'CLRSCRN’ ) 


BEAM ENERGY SELECTION: USED TO CALCULATE BETA & GAMMA 


oe le) 


PRINT *, ‘SELECT BEAM ENERGY: ' 
PRINT *,'1=50 MEV 5 2= SELECT YOUR OWN’ 
READ *,IBEAM 
IFC IBEAM .EQ. 1)THEN 
EBEAM=50. 
GO TO 300 
ELSE 
PRINT *,°INPUT THE DESTRED ELECTRON BEAM ENERGY IN MEV' 
PRINT *,°C INCLUDE DECIMAL POINT )' 
READ *,EBEAM 
END IF 
300 CONTINUE 
GAMMA=EBE AM/EREST TRIO0100 
BETA=SQRT(1.-( 1. /GAMMA#2. )) TRIOO110 
CALL EXCMS ( ‘'CLRSCRN’ } 


PRINT *, "SELECT A FREQUENCY BAND TO BE SUMMED OVER AND THE’ 
PRINT *,‘INCREMENT IN MHZ: 1 = A SUM OF FREQUENCIES' 
PRINT *,'FROM 10 MHZ TO 100 MHZ IN 10 MHZ INCREMENTS‘ 
PRINT *,‘'2 = SELECT YOUR OWN’ 
READ *, INCREM 
IF (INCREM .EQ@. 1)THEN 
BNU=10 .E6 
DLOOP=9 
DNU=10.E6 
CALL EXCMS ('CLRSCRN' ) 
GO TO 400 
ELSE 
CALL EXCMS ( 'CLRSCRN’ ) 
eccecceccceccceceecceccececcceececceeccceccecccecceccceccecceccccccccccccccccccce 
Cc CALL THE SUBROUTINE USED TO DEFINE FREQUENCY BAND OF INTEREST C 
Cc AND THE INCREMENTAL STEPS FOR SUMMING Cc 
eccececcececcecececcecececcceccccceccccccecccceccccecccceccececccecccceccceccecccccccc 
CALL FREQ( BNU,DNU,DLOOP ) 
END IF 
400 CONTINUE 
c 
Cc OUTPUT GRAPH AXIS SELECTION 
PRINT *,'SELECT OUTPUT TYPE’ 
PRINT *,° L=ENERGY VS. FREQUENCY FOR A SELECTED ANGLE THETA’ 
PRINT *,'2=ENERGY( INTEGRATED OVER FREQ''S OF INTEREST) VS. THETA’ 
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PRINT *,'3=FORMFACTOR(F<K>) VS. FREQ. FOR A SELECTED ANGLE THETA' 
READ *,SELECT : 
CALL EXCMS ( ‘CLRSCRN' ) . 
IF (SELECT .EQ. 1 .OR. SELECT .EQ. 3) THEN 
GO TO 470 
ELSE 
GO TO 450 
END IF 
450 CONTINUE 
PRINT *, ‘SELECT A RANGE FOR THETA AND THE INCREMENT AT WHICH ' 
PRINT *, 'THE RADIATED ENERGY WILL BE CALCULATED' 
PRINT *, ‘1 = 0 TO 180 DEGREES IN ONE DEGREE INCREMENTS' 
PRINT *, '2 = SELECT YOUR OWN' 
READ *, NUMBER 
IF (NUMBER .EQ. 1)THEN 


D2LOOP=181 
DTHETA=0.0 
DEGINC=1 
ELSE 
cccceccccccceccccccecccccceccceccecccccceccccecccccccccccccccccccccccccccccccc 
C CALL THE SUBROUTINE USED TO DEFINE THE INCREMENTAL C 
C STEPS OF THETA C 


ccccccecceccccccccccecceccecececcccccccccccce cece ccc ccccccccccccccccccccccccc 
CALL DEGREE (DTHETA,DEGINC ,D2LOOP ) 
END IF 
G65 CONTINUE 
GO TO 500 
470 PRINT *,'SELECT AN ANGLE (THETA) OF INTEREST. ' 
READ *,DTHETA 
MIDNU = BNU 
TTERG = 0.0 
DLOOP=DLOOP +1 
GO TO 510 
500 CONTINUE 
DO 600 K=1,D2LO0P TRIO0140 


C THIS OUTER LOOP WILL BE USED ONLY FOR OUTPUT TYPE 2 
TTERG=0. 
MIDNU= BNU+(DNU/2. ) 
510 CONTINUE 
THETA=DTHETA*PI/180. 
COSANG = COS( THETA) 
KONST = 2.%SIN( THETA )/((1./( IND*BETA ) )-COSANG) 
ATHETA = ( LENGTH¥PI/C )%((1./( IND¥BETA ) )-COSANG ) 


C THIS INNER LOOP IS USED FOR ALL 3 OUTPUT TYPES BY INCREMENTING 
C FREQUENCY OVER THE CHOSEN RANGE 
DO 550 J=1,DLOOP 
IFC IFORM .EQ. L)THEN 
GO TO 515 
ELSE 
COSARG 2.*BEE*PI/C 
COSGAR 2. *DEEXPI/C 
KKONST (C¥*2. 7102.40 PLZ. DRC (I DEEREZ. )-(BEEM*Z. ))) 
FORFAC=KKONST (1. /(MIDNU%#2. ) 301.70 ABSCCOSANG ) 2. ) J 
D (COS( COSARG*COSANGXMIDNU )-COS( COSGAR*COSANG#MIDNU ) ) 
END IF 
515 CONTINUE 
RADFUN=ABS( KONST *SIN( ATHETA#MIDNU )*¥FORFAC ) TRIOO180 
TERG=CONST( RADFUN##2. ) 
IF (SELECT .EQ. 1) THEN 
WRITE( 21,520) MIDNU,TERG 
520 FORMAT(2E10.5) 
ELSE IF (SELECT .EQ. 3)THEN 
WRITE (21,525) MIDNU,FORFAC 
525 FORMAT(E10.5,E10.¢) 
ELSE 
TTERG = TTERG + TERG 


530 CONTINUE 
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END IF 
MIDNU= MIDNU+DNU 
550 CONTINUE TRIOOZOO 
IF(SELECT .EQ. 1 .OR. SELECT .EQ. 33 THEN 
GO TO 600 
ELSE 
GO TO 555 
END IF 
555 CONTINUE 
ERG(K J=TTERG® DNU 
GEE =ABS( KONST ) #2. 
WRITE( 21,560) DTHETA,ERG( K } TRIOO190 
560 FORMAT (2E10.5) 
DTHETA = DTHETA + DEGINC 
600 CONTINUE 
C THE DATEND STATEMENT IS USED TO END THE DATA FILE 
Cc (NECESSARY IF USING DISSPLA GRAPHICS) 
DATEND=999.99 
WRITE( 21,660) DATEND,DATEND 
660 FORMAT (2F10.2) 
PRINT *, ‘DO YOU WANT TO MAKE ANOTHER RUN?(Y/N}° 
READ (5,670 )RESPON 
670 FORMAT(A1 ) 
IF(RESPON.EQ. 'Y'} THEN 
CALL EXCMS ('CLRSCRN’ ) 


GO TO 100 
ELSE 
GO TO 700 
END IF 
700 CONTINUE 
STOP TRIOO2Z10 
END TRIOO220 
Cc 
Cc 


cecccceccecccccecceccccececcccceccccccecccccecccccccccccccccccccccccccccccece 
Cc 
SUBROUTINE FREQ(BFREQ,DNU,DLOOP } 

REAL BFREQ,EFREQ,DIFF ,DNU 

INTEGER DLOOP 

PRINT® ,'ENTER BEGINNING FREQUENCY IN MHZ( INCLUDE DECIMAL POINT }' 

READ*® ,BFREQ 

BFREQ=BFREQ¥1.E6 

PRINT® ,'ENTER ENDING FREQUENCY IN MHZ( INCLUDE DECIMAL POINT }' 

READ® ,EFREQ 

EFREQ=EFREQ¥1.E6 

PRINT® ,°ENTER DESIRED INCREMENT IN MHZ( INCLUDE DECIMAL POINT }' 

READ® ,DNU 

DNU=DNU*1.E6 

DIFF=EFREQ-BFREQ 

DLOOP=INT(DIFF/DNU} 

RETURN 

END 
Cc 
Cc 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C 

SUBROUTINE DEGREE( BDEG,DEGINC ,DZ2LOOP } 

REAL BDEG,EDEG,DEGINC ,DIFFF 

INTEGER D2LOOP 

CALL EXCMS (‘'CLRSCRN' } 

PRINT *,'ENTER STARTING AND ENDING ANGLES OF INTERESTCREAL})' 

READ *,BDEG,EDEG 

PRINT *, ‘ENTER THE INCREMENTAL VALUE OF THETA( DEGREES) DESIRED. ' 

READ *, DEGINC 

DIFFF = EDEG-BDEG 

D2LOOP =INT(DIFFF/DEGINCG J+41 

RETURN 

END 
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